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Abstract 

The secret cave cricket, Ceuthophilus secretus Scudder (Orthroptera: Rhaphidophoridae), is an 
obligate trogloxene endemic to central Texas, USA, and is a primary source of energy and nutri- 
ents for sensitive cave ecosystems. In this study, nine polymorphic microsatellite markers were 
developed from genomic DNA of C. secretus. Genotypes of 120 individuals sampled from four 
localities on the Fort Hood Military Reservation, Killeen, Texas, were analyzed to characterize 
the polymorphism at each locus. The number of alleles ranged from 19 to 46. All paired loci were 
in linkage equilibrium. Four loci had significant deviations from Hardy- Weinberg equilibrium. 
Observed heterozygosity varied from 0.793 to 0.959. These loci provide a means of characteriz- 
ing population genetic structure in this species. 
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Introduction 



The Edwards Plateau of central Texas is one 
of the largest continuous karst regions in the 
USA, subsidence sinkholes and caves being 
the major karst features of this area (Venni 
1994). The many small caves located 
throughout this region contain unique ecosys- 
tems, with some harboring a variety of rare 
endemics, including several endangered inver- 
tebrates (U.S. Fish and Wildlife Service 1988, 
1993, 2000). Many of these species are tro- 
globites that do not leave the cave 
environment. The secret cave cricket, 
Ceuthophilus secretus Scudder (Orthroptera: 
Rhaphidophoridae), which is found in many 
cave ecosystems in central Texas, is consid- 
ered to be a trogloxene (Taylor et al. 2005; 
Lavoie et al. 2007) because it forages on the 
surface at night, but typically returns to the 
cave for other parts of its daily activities. En- 
ergy and nutrients, in the form of C. secretus 
biomass, feces, and eggs, are thereby trans- 
ferred from the surface into the cave 
ecosystem, where they are utilized by various 
predators and scavengers and become part of 
cave nutrient and bioenergy cycles (Taylor et 
al. 2007). Given the species' importance to the 
cave ecosystems, it is considered a keystone 
species on which many of troglobites depend 
(Taylor et al. 2005; Lavoie et al. 2007). 

Genetic markers are largely lacking for the 
family Rhaphidophoridae, excepting five mi- 
crosatellite DNA loci described for use in the 
genus Dolichopoda (Bernardini and Ketmaier 
2002). Nine new microsatellite DNA loci that 
can be used in characterizing population ge- 
netic structure in C. secretus were developed 
in this study. These markers should assist re- 
searchers in characterizing gene flow, and 
general dispersal and colonization patterns, 
within and among C. secretus populations. 
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Landscape and connectivity (dispersal, migra- 
tion, gene flow) are likely important features 
of C. secretus natural history, as evidenced by 
the species' highly variable residency and 
abundance in different caves over time (C. 
Pekins, Environmental Division, US Army 
Garrison-Fort Hood personal communication). 
Therefore, the ability to use genetic data to 
elucidate connectivity among caves (or popu- 
lations of multiple caves) and to identify 
landscape features that influence connectivity 
should prove useful for the management and 
conservation of this critical player in cave 
ecology. 

Materials and Methods 

Whole genomic DNA was extracted from the 
legs of each C. secretus specimen using a 
DNeasy Blood and Tissue Kit following the 
manufacturer guidelines (Qiagen, 

www.qiagen.com) . By utilizing the DNA 
from eight individuals, Genetic Identification 
Services (GIS; www.genetic-id-services.com) 
was able to develop four enriched microsatel- 
lite genomic libraries. These libraries were 
enriched for CA, AAC, ATG, and TAGA mo- 
tifs. GIS designed polymerase chain reaction 
(PCR) primers for 77 microsatellite- 
containing clones using Designer PCR version 
1.03 (Research Genetics, Inc., 
www.lifetechnologies.com) . From this pool of 
potential loci, reaction conditions were opti- 
mized for nine loci showing evidence of 
polymorphism. 

Microsatellite amplifications were performed 
in 10 ul reaction volumes with 1.0 ul of DNA 
template in a solution of 0.18 mM of each 
primer, 1 X NH 4 PCR Reaction Buffer (Bio- 
line, www.bioline.com) , 0.80 mM dNTPs, 
0.025 units of Biolase DNA polymerase (Bio- 
line), and 0.25 mM MgCl 2 . All PCRs were 
performed in a MyCyler Thermal Cycler (Bio- 
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Rad, www.bio-rad.com) under the following 
conditions: initial incubation at 94° C for 3 
minutes, followed by 30 cycles of 94° C for 
40 seconds, 54° C to 59° C (according to each 
primer annealing temperature; see Table 1) for 
40 seconds, and 72° C for 30 seconds, fol- 
lowed by a final 10 minute incubation at 72° 
C. 

In order to characterize microsatellite loci, 30 
individuals from each of four populations (N 
= 120) located in the Fort Hood Military Res- 
ervation near Killeen, TX, USA, were 
genotyped. The four populations, from Geo- 
cache Cave, Lost Chasm Cave, New Cave, 
and Sanford Pit Cave, were located on four 
different limestone ridges and separated by an 
average distance of 18.8 km (SD = 11.3 km). 
Individuals were genotyped using an ABI 
3130 Genetic Analyzer with labeled forward 
primers (6-FAM, NED, HEX) and ROX 400 
size standard (Applied Biosystems, 
www.appliedbiosystems.com) . Each locus 
was run singly on the genetic analyzer. All 
loci were tested for linkage disequilibrium and 
deviations from Hardy- Weinberg expectations 
using GENEPOP 4.1.0 (Rousset 2008). A se- 
quential Bonferroni correction (Rice 1989) 
was applied to correct type I error rates for 
multiple significance tests. Observed and ex- 
pected heterozygosities were calculated by 
using GENALEX 6.0 (Peakall and Smouse 
2006). GENEPOP was used to test for signifi- 
cant differentiation in allele frequencies 
among populations and to calculate interpopu- 
lation FSTs. 

Results 

Genotyping results are summarized by locus 
in Table 1. All nine loci were polymorphic, 
with numbers of alleles ranging from 19 to 46. 
Observed heterozygosity ranged from 0.793 to 
0.959, while expected heterozygosity values 



ranged from 0.854 to 0.969. Differences be- 
tween observed and expected heterozygosity 
were due, at least in part, to having estimated 
heterozygosity using data pooled from four 
genetically differentiated populations (the 
Wahlund effect). Four of nine loci showed 
significant deviations from Hardy-Weinberg 
equilibrium (HWE) in at least one population. 
MICROCHECKER 2.2.3 (Van Oosterhout et 
al. 2006) analysis suggested that the presence 
of null alleles might explain heterozygote de- 
ficiencies for Cese2 (H E = 0.957, H 0 = 0.704, 
null allele frequency (NuIIf) = 0.132) and 
Cese9 (H E = 0.910, H 0 = 0.704 Null F = 0.108) 
in POP1, for loci Cesel (H E = 0.952, H 0 = 
0.844, NuIIf = 0.057) and Cese3 (H E = 0.931, 
Ho = 0.719, Null F = 0.115) in POP2, and for 
locus Cese2 (H E = 0.945, H 0 = 0.793, Null F = 
0.081) in POP4. There were no cases of link- 
age disequilibrium among loci. All allele 
frequencies were significantly different be- 
tween pairs of populations (p < 0.009), and 
interpopulation FSTs ranged from 0.001 to 
0.014. 

Discussion 

In this study, a suite of microsatellite DNA 
markers were developed for the secret cave 
cricket, C. secretus. These markers should 
provide insight into genetic diversity, gene 
flow, and landscape genetics in a species for 
which little such information is available. In a 
species such as C. secretus, which resides in a 
habitat that is difficult to access, with highly 
unpredictable presence and abundance among 
caves and over time, and that is largely noc- 
turnal and secretive (hiding under rocks and 
debris) when on the surface, molecular tech- 
niques likely provide the best means for 
determining colonization and dispersal pat- 
terns. 
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Four of nine microsatellite markers showed 
genotypic frequencies that deviated from those 
frequencies expected under HWE. Population 
substructure, inbreeding, or technical artifacts 
(e.g., null alleles) could be the cause for these 
deviations. Additional populations should be 
examined to determine if these loci consistent- 
ly show evidence of null alleles and deviations 
from HWE. No evidence of deviations from 
HWE expectations for the other five loci was 
found. Although care should be used when 
including loci with null alleles in some anal- 
yses, the resulting bias is often small 
(Carlsson 2008), and methods exist allowing 
their inclusion in several estimators of genetic 
parameters (van Oosterhout et al. 2004; Cha- 
puis and Estoup 2007). 
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Table I. Locus name followed by repeat motif present in sequenced fragment, primer sequences with designated fluorescent dye, 
primer annealing temperature (Ta), and size range of alleles in base pairs (range). Also, number of alleles (N a ), observed heterozy- 
gosity (Ho), and expected heterozygosity (He) based on combined data from four populations of Ceuthophilus secretus are reported. 
Differences in Ho and He can be at least partly explained by the Wahlund Effect. 



Locus 


Repeat Motif 


Primer Sequence 


T A 


Ranee 


N A 


Ho 




Cesel 


(CA) 20 


F(6FAM): CCTTGACAGTTGAACTCAGTG 


55°C 


122-250 


46 


0.861 


0.969 


R-TTTCTTGAACGGGATATACAA 


Cese2 


(CA)„ 


F(HEX): TTCCATGCGAAATAGATGACC 


55"C 


223-303 


37 


0.814 


0.963 


R-CGAACCAGCGTACTACGAGTC 


Cese3 


(CA), 4 


F(NED): GCACGCAGTAGAGATTATCC 


59° C 


284-352 


33 


0.831 


0.914 


R-CCCCTGTAACAAGAAATCATT 


Cese4 


(AAC),, 


F(6FAM): GGCTACTGTTATTGCTTCTCG 


54° C 


227-290 


19 


0.833 


0.961 


R-GCTTAGGTTGCGAACACTAAC 


Cese5 


(AAC) 14 


F(NED): CTCCGCTCGGTGTAACAG 


54° C 


250-334 


21 


0.826 


0.920 


R-CATACGCAGTAAGCCAAAAAT 


Cese6 


(ATG) 7 


F(HEX): CACAACTGAGGCGATTAAGAC 


57° C 


248-320 


20 


0.902 


0.901 


R-GAACGGGACTCAAAGGTTATC 


Cese7 


(TC)„ 


F(NED): TCTGAATCGCCTCAGTATGA 


55" C 


115-207 


35 


0.959 


0.857 


R-CACACGAAGGAGGTTGAAG 


Cese8 


(TAGA) 25 


F(6FAM): GAAGGCAGCAACTTTTCAC 


59° C 


100-348 


41 


0.943 


0.854 


R-AGGGTGGCAATGTTTAGG 


Cese9 


(TAGA) 6 


F(NED): ACATCAACTTCCTACGAGTGTC 


54° C 


124-268 


26 


0.793 


0.911 


R-AAGATTTTGGCTTCAGTCAAG 
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